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Abstract. We study a lattice gauge theory in Wilson’s Hamiltonian formalism. In 
view of the realization of a quantum simulator for QED in one dimension, we introduce 
an Abelian model with a discrete gauge symmetry Z„, approximating the 17(1) theory 
for large n. We analyze the role of the finiteness of the gauge fields and the properties 
of physical states, that satisfy a generalized Gauss’s law. We finally discuss a possible 
implementation strategy, that involves an effective dynamics in physical space. 
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1. Introduction 

Recent development in low-temperature physics and atomic control techniques is 
providing the basic tools for setting up quantum simulators The experimental 

feasibility of a quantum simulator will open the way to a more comprehensive 
understanding of complex systems and fundamental physics. An appealing application 
is the simulation of lattice gauge theories [a El m El El [In] . The discretization of 
high-energy physics theories on lattices was initially motivated by the possibility to 
simulate them by classical computation. However, the complex nature of gauge theories 
represents a severe obstruction, which can be overcome through a quantum simulation. 
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Atoms on a lattice [HI [T2l [I3l [HI [151 CBl IHl [HI [IHl EQl ET] provide a natural toolbox 
to perform this task. 

In this article we construct a lattice model which simulates an Abelian gauge 
theory. The model will be characterized by the interaction between an Abelian gauge 
(electromagnetic) field and a fermionic matter field, as in quantum electrodynamics 
(QED). We shall restrict our analysis to the one-dimensional case. As in the consolidated 
quantum link model (QLM) [221 ESI EH EHi system will be an approximation to 
QED in which the electric field can take a finite number of values. We will eventually 
discuss the possibility to implement the model on a cold atomic simulator. This task 
involves first the identification of the degrees of freedom of the simulator with those of 
the model, and then the correct implementation of the dynamics. 

2. Lattice QED 

The dynamics of the continuum-space QED in 1 -|- 1 dimensions (Schwinger model) has 
its lattice counterpart in the Hamiltonian [6l ESj EE] 

X XX 

with X labelling the sites of a one-dimensional lattice. Fermionic matter is represented 
by the one-component spinor field operators fjx, defined on each site, which obey the 
canonical anticommutation relations {V’*,'*/'!'} = ^x,x', = 0. The parameter 

m is the fermion mass, while the staggered fermions with parity factor (—1)* are 
introduced in order to avoid the fermion-doubling problem in the discretization of the 
theory p El El [E|: the positive and negative-energy components of the Dirac spinor 
are encoded respectively in the even and odd lattice sites. In a simplified model, 
spinless particles are considered, with the spinor reducing to a single-component 
held. The gauge helds are instead dehned on the links (x, a; -1- 1) of the lattice. In 
the canonical gauge, which is the most convenient choice to develop a lattice gauge 
theory, the electric held E and the vector potential A are conjugated variables with 
canonical Commutation Relations (CR) [Ex^x+i, Ax\x'+i\ = ihx,x' [22|- Since the unitary 
operators U (comparators [2H]) in the hopping terms of ([^ are locally related to the 
vector potential by exponentiation, Ux,x+i = ^ they satisfy 

[h/XjX+i) Ux',x'+i\ hx,x'Ux,x+i- (2) 

The presence of the free electric held energy, with coupling constant thus yields a 
nontrivial dynamics for the comparators. Notice the absence of magnetic contributions 
to the Hamiltonian, a consequence of the one-dimensional nature of the system. 

The terms 'ilj'lXJx,x+i'f’x+i in ([^ describe site hopping of fermions, related to a shift 
in the electric held. These contributions come from the discretization and integration 
on the lattice cells of the minimal-coupling terms —i' 0 (a;)i 7 o 7 -^Dj' 0 (x), with 7 ^ the 
Dirac matrices, and Dj = dj + iAj{x) the covariant derivatives P EH]. The minimal 
coupling ensures the symmetry of the Hamiltonian under local U{1) transformations. 


Discrete Abelian Gauge Theories for Quantum Simulations of QED 


3 


Given a real function on the lattice, ax, local phase transformations xjjx —)■ and 

Ux, x+i —>■ D^'^^Ux^x+i^ of the held operators leave the Hamiltonian Q invariant. 
Due to the (anti)commutation properties of the helds, the phase transformation on any 
operator F can be implemented through the application of 

F (3) 

X y 

Gx = 'ifl'lfx + ^[(- 1 )" - 1 ] - {Ex,x+1 - Ex-l,x). 

The gauge transformation acts trivially on the local-17(l) invariant Hamiltonian 
For a generic operator F, the invariance property is equivalent to the CR [F, = 0 

for all sites. The choice of the canonical gauge has the disadvantage that Gauss’s law 
cannot be implemented at the operator level [2ZII2S]. Thus, not all the states in the 
total Hilbert space TL of the system are physically acceptable. The invariance condition 
of a state |0) under any gauge transformation ~ 1*?^) selects the Hilbert 

subspace 

TLg = {|0) £ "E, Gx\(l)) = 0 for all sites x} , (4) 

where, as can be deduced from (|^, the charge V’iV’x and the divergence of the electric 
held Ex^x+i — Ex-i^x are correctly related. The term multiple of the identity appearing in 
Gx ensures that the vacuum state with vanishing electric held and all the negative-mass 
sites occupied (the Dirae sea) is in the gauge-invariant subspace TLg- 

3. Finite link spaces: a model. 

In Wilson’s original formulation of lattice gauge theories [7j, the gauge operators on 
links act on inhnite-dimensional Hilbert spaces, and both the electric held and the 
vector potential have continuum and unbounded spectra. A problem arises in quantum 
simulators, when one has to match the (inhnite dimensional) link with an experimentally 
feasible and controllable system, with a finite number of levels. Two possible approaches 
are possible. The hrst one preserves for all dimensions the structure of the Hamiltonian 
([^, including the coupling of the matter helds with a unitary gauge operator. The 
second approach, which has been followed in the formulation of the QLM [22j, consists 
in preserving the GR between held operators, (|^ in particular. Since in quantum 
mechanics one is used to think in terms of commutators, the latter approach has so far 
appeared more natural. Unfortunately, this procedure focuses on the invariance with 
respect to the U(l) group of local transformations, at the expenses of the structure 
of the hopping term, that no longer involves a unitary comparator (minimal coupling 
prescription). We shall rather insist on the unitarity of the comparator thereby obtaining 
a bona fide lattice gauge theory for any dimension of the link Hilbert space. Let us hrst 
observe that a gauge transformation ([^ acts on the comparator Ux ,x+i = e as 

TT _\ — 0'x+l)Ex,x+l T T ^ — (A'x-\-l)Ex,x-\-l —Q^cc+l)/’/’ 


( 5 ) 
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This result is indeed a special case of a general property of operators in the Weyl group 
generated by the conjugated operators A and on a link (indices will be omitted for 
clarity) [2911301 • Indeed, the electric field and the vector potential are the generators of 
the two-parameter projective unitary (Weyl) group Using the canonical 

CR [E, A] = i for the generators, the following relation holds for any ^ and r] 

which particularizes to ([^ for r; = 1 and f 

In a hnite-dimensional Hilbert space, the role of generators loses its meaning. 
Nonetheless, we can define a set of unitary operators, the (discrete) Schwinger-Weyl 
group, that satisfy the relation ([^. Notice that this entails a change of paradigm. 
We are abandoning an approach in terms of the algebra of generators in favor of an 
alternative approach in terms of its group, enabling us to explore the global features of 
the topology. This concept is elaborated in the Appendix. Observe that the operator 
e-ir)A ^ translation of the electric field, since = E + rj. A similar role 

is played by the elements of the Weyl group with ^ = 0 with respect to A. 

Let us now consider an n-dimensional Hilbert space and choose an orthonormal 
basis {|nfc)}i<fc<,fi, which will be called the electric field basis, and dehne a unitary 
operator U that performs a cyclic permutation of the basis states: 

U\vk) = \vk+i) for k <n, U\vn) = |ni). (7) 

We will call the orthonormal eigenbasis of U the vector potential basis. The operator V 
conjugated to U is diagonal in the electric held basis, with V\vk) = e~^'^^^^'^\vk), and it 
cyclically permutes the elements of the vector potential basis. Since = U” = I, the 
sets {U^}i<fc<n and {U^}i<fc<n are unitary representations of the group of integers 
modulo n. The set of all the products between U and V and their integer power 
constitutes the Schwinger-Weyl group [2S1 EDI El]- The CR between the elements of this 
group yield the relation 

V-ku^V^ = with k,lez. (8) 

The multiplication law (|^ satished by the operators is the discrete form of (|^, valid 
for U{1) operators. 

Once the correspondences ^ and V^^x+i ^ has been set 

up for all links, we can construct an Abelian theory which represents an approximation 
to the lattice QED Hamiltonian ([^ with a local invariance |32l ES]. One of the 
earliest examples of a pure gauge model with a invariance, used to approximate a 
U{1) theory in the n —)■ oo limit, was given in [M]. The dynamics of the new model is 
determined by the following variant of the Hamiltonian (|^ 

X XX 

where f{V) is a suitable Hermitian operator, diagonal in the electric field basis, which 
represents the discretized free-held electromagnetic Hamiltonian. Unlike in the QLM, 
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the ladder operators U permute the electric held basis states on a circle, and transitions 
between neighboring states all occur with the same amplitude. These features are a 
consequence of the request that the minimal-coupling structure in Q is preserved in 
its hnite-dimensional link counterpart. The function f{V) is so far arbitrary. We shall 
introduce a function f{V) that, like E'^ in Q, has a single minimum. A simple choice 
is 


f(y^x,x+i) ^{^x,x+i 


i)(n',.+i -1)- 


( 10 ) 


the operator (10) has eigenvalues 


Since the eigenvalues of V are Vk = 

S(kx,x+i) = (sin low-energy states around the minimum at fc = 0, the 

spectrum of f(V) is quadratic like the energy E^ associated to the electric held in the 
original model. See Fig. [Tj We remark that the choice of a function for the electric held 
energy does not ahect the discretization of the gauge theory, as long as it is characterized 


by a single minimum. The energy term (10) is characterized by a nondegenerate ground 
state. Another interesting possibility is obtained by the replacement V —)■ ^ 

yielding two degenerate minima at fc = — 1 and /c = 0, which could be of interest for 
studying the ehects of terms that break the chiral symmetry of the gauge held |35l |36l Il2] . 
As expected, in the link discretization process the U{1) gauge symmetry of the original 
theory becomes a symmetry. Indeed, due to condition (|^, the Hamiltonian Hn is 
invariant under the transformation 


Hn —t 






( 11 ) 


where Vx is an arbitrary integer-valued function on the lattice, and 


T 

-L (7 




2 


'^X,X+lVx-l X-: 


( 12 ) 


with (Ta;)" = I, represents the discretized analogue in the Schwinger-Weyl theory of 
[see Eq. ([^], obtained by the correspondence e“‘® —)■ V. Due to the arbitrariness of z/j., 
the gauge-invariance condition is equivalent to 

lHn,Tx]=0, Vx. (13) 

In the Zn theory, the Hilbert subspace of physical states is determined by a generalized 
Gauss’s law: 

Ht = {\4>) £ H, Tx\(j)) = 10) for all sites x} . (14) 

In the reference basis {\nx, W,x+i)}, determined by the eigenvalues of 14,a;+i and of 
namely the occupation numbers Ux € {0,1}, the gauge invariance request translates into 
a condition on the eigenvalues. Since each Tx acts nontrivially on site x and the adjacent 


links, the physical subspace (14) is spanned by the reference basis states satisfying 


^(n 
e n V 


(-i)^-i j 

’ = 1. Vi. 


( 16 ) 


If an even site with (—1)’” = 1 is empty, the eigenvalues of V in neighboring links must 
be equal, while if it is occupied, they will be related by W,x+i = e On 
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Figure 1. Spectrum S{kx,x+i) of the operator f{Vx^x+i) defined in (10). The dashed 
(red) line represents the quadratic approximation of the spectrum close to its minimum 
at fc = 0. 


the other hand, in odd sites with (— 1 )"^ = — 1 , the eigenvalnes of V are eqnal if the 
site is occnpied, while they are related by Vx,x+1 = e^^^/^Vx-l,x otherwise. This is in 
agreement with Dirac’s picture in which the absence of a particle in a negative-energy 
site is equivalent to the presence of an antiparticle. We can represent these situations 
by visualizing the eigenstates of each V as n points placed at an angular distance of 
2Ti/n on a circle. Figure [^displays the state of an occupied even site. The action of the 
operator Ux,x+i (f^x+i) induces a counterclockwise (clockwise) jump of the eigenstate 
of D to a neighboring state on the circle. The correlation between the jump in the 
link eigenstate and the charge displacement in the terms ipl.Ux^x+i'f’x+i (and conjugates) 
keeps the system in the physical subspace TLt- 

Although the approximation of Wilson’s model with a lattice gauge theory 
improves with increasing n, interesting phenomenology emerges also at small link 
dimension: the Z 3 theory is already a good test bed for the analysis of electric flux 
string breaking [ISIET]. Moreover, theories are relevant per se, being related to the 
problem of conhnement in QCD |33] . 

Let us hnally remark that an extension to higher dimensional lattices would involve 
a plaquette term in the Hamiltonian with a product of four U matrices [25]. The 
theory can be generalized by a proper modihcation of Gauss’s law, which takes into 
account all the terms in the divergence of the electric held. 
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Vx,x+\ 


Figure 2. The state of two links (large circles) neighboring an even site [central full 
dot, p = (—1)^ = !]■ In this example, the site is occupied. A full dot on the circle 
represents the occupied eigenstate of V, while the empty dots are the other eigenstates. 
Observe that unitarity entails probability conservation: the full dot never leaves the 
circle under 0. 

4. Implementation strategy 

The implementation of the Hamiltonian ([^ is a difficult task due to the presence of the 
correlated hopping terms describing elementary processes in which the 

hopping a: —)■ x + 1 of a fermion to a nearest-neighboring site is always associated with an 
electric-field change on the link {x,x + 1) between the two sites. Recent results indicate 
that it is possible to engineer a simpler Hamiltonian, and then obtain the gauge theory 
by imposing the gauge-invariance constraint (s) by assigning an energy penalty to the non 
gauge-invariant subspace [121 EH] or by implementing Zeno constraints [211EU]EUIITTI E2] . 
We shall follow the first path, and induce the correlated hopping in (|^ through a second- 
order effective gauge-invariant Hamiltonian. 

In view of implementing a Hermitian term in the Hamiltonian, it is more convenient 
to express the local gauge-invariance condition = |0) as 

r.|0)=o, r. = (T.-i)(T]-i), (16) 

so that the physical space is the kernel of the positive operator 

r = (17) 

X 

The dynamical conservation of the gauge condition for all sites can be enforced by 
adding a large term proportional to T to the Hamiltonian. Notice that the constraint is 
diagonal in the reference basis {\nx,Vx^x+i)}- 

Let us consider a Hamiltonian that involves uncorrelated hopping of fermions 
between nearest-neighbor sites and transitions, with equal amplitudes, between neighbor 
link states on a circle: 

^(0) = + Hd. (18) 
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The part Hd includes all the terms that are diagonal in the reference basis, such as the 
fermion mass term the gauge field energy (fi'n/ 2 ) /(14_3;+i), and 

possibly proper counterterms. 

The fermion hopping terms in (18) emerge naturally in condensed matter physics 
when a tight-binding approximation is assumed, in which only the lowest band of the 
lattice is energetically accessible. The matter field can be represented by a fermionic 
atomic species: the operator ‘ipl creates an atom in the fundamental Wannier function 
centered on the lattice site x [13]. Tunneling between neighboring sites provides the 
mechanism for the hopping processes The staggered structure is obtained by 

modulating the depth of the lattice wells [12] . A physical implementation of the action 
of the operators U on links is less obvious, since it is first of all necessary to identify 
a proper system in which transitions between adjacent levels on a circle occur with 
the same amplitude (see Fig. |^. A possible implementation tool is represented by a 
longitudinal array of transverse ring-shaped lattices [H], with their axes aligned with 
x, each one representing a link and conhning a single boson or fermion (the statistics 
being immaterial, since each ring contains one particle). The system is represented in 
Fig. In the tight-binding regime, the fundamental Wannier states, centered on ring 
lattice sites, can be identified with the link reference basis states |na;,a;+i)- Link particles 
are associated with the field operator C(^i-)x,x+i, with i = 1,... ,n the site index. Hopping 
between neighboring sites on the link yields the processes 


a 


x,x-\-l 




2=1 


2=1 


with C(„+i) = C(„). Transition between different rings are instead forbidden by a large 
energy barrier. Thus, the system and the hopping mechanism provide both the circular 
structure and the equality of transition amplitudes. Other possible implementations 
could involve a coupling between the internal levels of an atom confined on the link. 
However, in this case the unitary properties of transitions do not emerge by symmetry, 
and would require a fine tuning of the transition amplitudes |45j . 

While the term Hd does not couple orthogonal eigenspaces of F^,, each nondiagonal 
term in (18) maps a state |0) of the physical space Ht = kerF into another eigenstate 
10') of each F^, such that 

F|0') = 27,10'), (20) 

where 7, = 2(1 — cos(27r/n)) is the first excited eigenvalue of F^,, and the factor 2 
is related to the fact that each fermion hopping or Rabi transition on a link affects 
the eigenvalue of F^. in two neighboring sites. The implementation strategy consists in 
adding to iL(o) a term which induces a large energy cost for going out of the physical 
space, thus obtaining 

H{i) = +'jn^uT. ( 21 ) 

If u is much larger than the parameters appearing in Ff{o), the evolution of an initial 
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Figure 3. Scheme of the physical implementation of the model. The fermions are 
confined in longitudinal lattice sites [light (red) spots], while the links are represented 
by ring-shaped transverse lattices, each one hosting one boson or fermion. The link 
particle is tightly confined in the longitudinal direction, while it can hop between 
nearest-neighbor sites [dark (blue) spots] on each ring lattice. 


state in the physical subspace TLt is approximately given by the effective Hamiltonian 


i/eff = PHa 


po)P - Pi/(0)Q(7n '«r)-^Qhf(0)P, (22) 

with P the projection operator on TLt and Q = I — P. Using the fact that PH{q)P = 
PHdP, and observing that P(o) couples TLt only with the eigenspace of V belonging to 
27 „, Eq. (22) simplifies to 

H.g = PHiP - (2u)-THf„fQHf„,P 


(23) 


A straightforward computation of the second term yields, up to immaterial constants, 
the final result 


Peff = P 


Pd- —+ H.C.)+— 


n„+i)jp, (24) 

in which the correlated hopping terms in (|^ appear. The counterterms in can be 
tuned to cancel the undesired (though gauge-invariant) contribution in (24). Despite 
the necessity to couple each site with the neighboring links through Tj., this strategy can 
be more convenient than a direct implementation of ([^, since the couplings between 
sites and links are all diagonal in the reference basis. It should be also noted that the 
operator T can be replaced by any positive operator whose kernel is TLt- 


5. Conclusions and outlook 

We have discussed a lattice gauge model with a discrete Abelian symmetry. The 
model represents an approximation to Wilson’s lattice QED, that improves as the 
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dimension n of the link Hilbert space increases. The obstacles to an experimental 
implementation of the model can be overcome through the use of effective dynamics. 
Further research will be devoted to two main avenues. The hrst one is the identihcation 
of an atomic or condensed-matter system which is suitable for an experimental 
implementation. The second is the application of the techniques developed in this 
Letter to non-Abelian gauge theories. Moreover, we shall scrutinize the convergence of 
the models towards U{1), also in comparison with the QLM, and analyze in detail 
the phenomenology of the model, including chiral symmetry breaking and conhnement. 
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Appendix A. Weyl group representation in infinite dimensions 


(A.l) 


Gonsider two Hermitian operators A and E with eigenvalue equations 

E\e) = e|e), A\a) = a\a), 

and satisfying the canonical commutation relations 

[E, A] = il. (A.2) 

Observe that the above commutation relations make sense in an inhnite-dimensional 
space. Two operators satisfying Heisenberg’s commutation relations generate a 
representation of the Weyl group. We are primarily interested in the unitary operators 

U(ri) := e-*’"’, V(0 := (A.3) 


with 77 G M, which act on the eigenstates (A.l) as spectral translations (in opposite 
directions): 

U{ri)\e) = |e + //), V(0|tt) = |a - 0- (A.4) 


Using the exponential form of the unitary operators (A.3), the commutation relations 


(A.l) and the Baker-Gampbell-Haussdorf formula in the form e^e'^ = e'^e^e^^’'^], valid 


if [R, S'] is a multiple of the identity, one obtains the property 

W(77)V(0 = e‘"«V(0W(h)- (A.5) 

The above result expresses the noncommutativity of E and A at the level of the Weyl 


group. The derivative of (A.5) with respect to f at (^, 77 ) = (0,1) yields 
[E,U]^U, 


(A.6) 


with U := W(l), which is Eq. (2) of the Letter. Notice that by taking the derivative 


of (A.5) with respect to r] and f at the origin one also reobtains Eq. (A.2). 
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In the infinite-dimensional case (A.5) implies (A.6). However, only (A.5) can be 


realized in the finite-dimensional case (for discrete values of rj and f) by preserving 


the unitarity of all the operators involved. By contrast, (A.6) does not admit a finite 


dimensional representation that preserves the properties of both E and U. Indeed, 
given two operators, U unitary and E Hermitian, on a finite-dimensional space, the 


commutation relations (A.6) cannot be satisfied, since they yield UEU"^ = E — I, which 


contrasts with UEU"^ being isospectral to E with a bounded spectrum. 

In the Quantum Link model of QED, which is based on the identification of the 
gauge degrees of freedom with spin variables E ^ Sz, U ^ S+ = + iSy, the unitary 

operator U is replaced with a non-unitary one in order to mantain the commutation 
relations [E, U] = U valid. We will instead choose to maintain the unitary structure 
of the gauge comparator U in the finite-dimensional case. To accomplish this task, we 
shall abandon the Heisenberg’s algebra relations between the operators A and E, which 
cannot be realized in finite dimensions. We will focus instead on the group relations 


between their complex exponentials (A.5), which admit a natural extension to finite¬ 


dimensional spaces through the representation of the Schwinger-Weyl group. 


Appendix B. The discrete Schwinger-Weyl group 

Let us consider an n-dimensional Hilbert space Tin and choose an orthonormal basis 
{\vi)}i<^<n- If is then possible to define a unitary operator U which rotates the basis 
states as 


U\vi) = \ve,+i) for £ < n, f/|n„) = |ni). (B.l) 

Since [/"■ = I, the integer powers of U constitute a representation of the group of 
integers modulo n. Note that the transition between the last and the first state of 


the basis [last line of (B.l)] is necessary to ensure the unitarity of U. The eigenvalue 
equation 

U\uQ=e-"-^^\uQ, 
with fc G {0,1 ,..., n}, is satisfied for 


(B.2) 


1 


W) = 

\/n ^ 


e " 




(B.3) 


r=i 


It is now possible to define an operator V which rotates the eigenbasis of E as 

V\uk) = \uk-i) for k> 1, V\ui) = \un). (B.4) 

The operator V is also unitary, with = I. From the definition (|B.4) and the form 


(B.3) of the eigenstates of U, one can easily demonstrate that {|'i^r)}i<£<„ is in fact the 


eigenbasis of V, with 

V\vi) = e~^^\ve). (B.5) 

The two operators V and U are thus called conjugated, since each one rotates the other 
one’s eigenbasis. It is also relevant to observe that the action of U and V on each 
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other’s eigenbasis in (B.1)-(B.4) is the hnite-dimensional counterpart of the spectral 
translations (A.4). (Notice the opposite signs.) The actions of U and 14 on a state do 
not commute. Indeed, comparing 

VU\uk)=e-"-^’^V\uk)=e-"-^^\uk-i) 

with 

UV\uk) = U\uk-i) = 

one obtains the relation 

UV = e^VU, 


(B.6) 


(B.7) 


(B.8) 


which can be immediately generalized to all the integer powers of the operators into 




(B.9) 


coinciding with Eq. (8) of the Letter. This result represents the finite-dimensional 
generalization of (A.5). The most striking difference is that the relation (B.9) is valid 
only for the discrete set of integers, which implies that the differentiation which leads 
from (A.5) to (A.6), is in this case meaningless. However, the procedure leading to 
(B.9) is successful in preserving the unitarity of the operators involved. Thus, (B.9) is a 
valid starting point for a gauge theory in which the local fields act on finite-dimensional 
Hilbert spaces. The continuum limit is recovered by introducing the Hermitian operators 


An and En, such that 


W = e 




yk _ g-gfeEn 


(B.IO) 


with r]i := l^jTnJn and := k^^2'n /n. The Hermitian operators satisfy (A.2) by taking 
the limit n —)■ cx), with rji ^ rj and fk ^ ^ [30] . 


References 

[1] Lewenstein M, Sanpera A and Ahufinger V 2012 Ultracold Atoms in Optical Lattices: Simulating 

Quantum Many-Body Systems (Oxford University Press, New York). 

[2] Cirac J I and Zoller P 2012 Nat. Phys. 8, 264. 

[3] Bloch I, Dalibard J and Nascimbene S 2012 Nat. Phys 6, 267. 

[4] Blatt R. and Roos C F 2012 Nat. Phys. 8, 277. 

[5] Rothe H J 1992 Lattice gauge theories (World Scientific). 

[6] Montvay I and Munster G 1994 Quantum Fields on a Lattice (Cambridge University Press, 

Cambridge). 

[7] Wilson K 1974 Phys. Rev. D 10, 2445. 

[8] Kogut J B and Susskind L 1975 Phys. Rev. D 11, 395. 

[9] Susskind L 1977 Phys. Rev. D 16, 3031. 

[10] Kogut J B 1979 Rev. Mod. Phys. 51, 659 (1979). 

[11] Kapit E and Mueller E 2011 Phys. Rev. A 83, 033625. 

[12] Banerjee D, Dalmonte M, Muller M, Rico E, Stebler P, Wiese U J and Zoller P 2012 Phys. Rev. 

Lett. 109, 175302. 

[13] Tagliacozzo L, Cell A, Orland P and Lewenstein M 2013 Nat. Commun. 4, 2615. 

[14] Zohar E, Cirac J I and Reznik B 2012 Phys. Rev. Lett. 109, 125302. 














Discrete Abelian Gauge Theories for Quantum Simulations of QED 


13 


[15] Kasamatsu K, Ichinose I and Matsui T 2013 Phys. Rev. Lett. Ill, 115303. 

[16] Banerjee D, Bogli M, Dalmonte M, Rico E, Stebler P, Wiese U J and Zoller P 2013 Phys. Rev. 

Lett. 110, 125303. 

[17] Tagliacozzo L, Celi A, Zamora A and Lewenstein M 2013 Ann. Phys. (Amsterdam) 330, 160. 

[18] Zohar E, Cirac J I and Reznik B 2013 Phys. Rev. Lett. 110, 125304. 

[19] Zohar E, Cirac J I and Reznik B 2013 Phys. Rev. A 88, 023617. 

[20] Hauke P, Marcos D, Dalmonte M and Zoller P 2013 Phys. Rev. X 3, 041018. 

[21] Stannigel K, Hauke P, Marcos D, Hafezi M, Diehl S, Dalmonte M and Zoller P 2014 Phys. Rev. 

Lett. 112, 120406. 

[22] Horn D 1981 Phys. Lett. lOOB, 149. 

[23] Orland P and Rohrlich D 1990 Nucl. Phys. B338, 647. 

[24] Chandrasekharan S and Wiese U J 1997 Nucl. Phys. B492, 455. 

[25] Wiese U J 2013 Annalen der Physik 525, 777. 

[26] K. Melnikov and M. Weinstein, Phys. Rev. D 62, 094504 (2000). 

[27] Fradkin E 2013 Field Theories of Condensed Matter Physics (Cambridge University Press, 

Cambridge, England). 

[28] Peskin M E and Schroeder D V 1995 An introduction to quantum field theory (Addison-Wesley 

Publishing Company, Reading). 

[29] Weyl H 1950 The theory of groups and quantum mechanics (Courier Dover Publications). 

[30] Schwinger J and Englert B G 2001 Quantum mechanics: symbolism of atomic measurements 

(Springer, Berlin). 

[31] Varadarajan V S 1995 Lett. Math. Phys. 34, 319. 

[32] Zohar E and Burrello M 2015 Phys. Rev. D 91, 054506. 

[33] Zohar E, Cirac J I and Reznik B 2015 Quantum Simulations of Lattice Gauge Theories using 

Ultracold Atoms in Optical Lattices arXiv: 1503.02312 [quant-ph]. 

[34] Elitzur S, Pearson R B and Shigemitsu J 1979 Phys. Rev. D 19, 3698. 

[35] Coleman S, Jackiw R and Susskind L 1975 Ann. Phys. 93, 267. 

[36] Coleman S 1976 Ann. Phys. 101 , 239. 

[37] Pichler T, Dalmonte M, Rico E, Zoller P and Montangero S 2015, arXiv:1505.04440 

[38] Zohar E and Reznik B 2011 Phys. Rev. Lett. 107, 275301. 

[39] Wu L A and Lidar D A 2002 Phys. Rev. Lett. 88, 207902. 

[40] Facchi P and Pascazio S (2008) J. Phys. A: Math. Theor. 41, 493001. 

[41] Wu L A, Kurizki G and Brumer P 2009 Phys. Rev. Lett. 102, 080405. 

[42] Raimond J M, Sayrin C, Gleyzes S, Dotsenko I, Brune M, Haroche S, Facchi P and Pascazio S 

2010 Phys. Rev. Lett. 105, 213601. 

[43] Bloch I, Dalibard J and Zwerger W 2008 Rev. Mod. Phys. 80, 885. 

[44] Amico L, Osterloh A and Gataliotti F 2005 Phys. Rev. Lett. 95, 063201. 

[45] Geli A, Massignan P, Ruseckas J, Goldman N, Spielman I B, Juzeliunas G, and Lewenstein M 

2012 Phys. Rev. Lett. 112, 043001. 




